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Abstract 

We derive the commutation relations for field strengthes of the non-Abelian gauge fields by re- 
quiring that the quantum mechanical equations of motion coincide with the classical field equations. 
The equations of motion with respect to time derivative coincide with the corresponding field equa- 
tions, while those with respect to space derivatives agree with classical field equations if constraint 
equations are fulfilled. Using the uncertainty relations for field strengthes at the same times, which 
follows from the commutation relations, we estimate the energy gap of vacuum fluctuations in the 
long wave limit as £o,fc ~ chhr 2 V~ l [g 2 /cK\ 1 ^ and the total energy £q ~ Ach (g 2 j 'hcV ^ with A 
being a cutoff separating the weak and strong coupling regimes. 

PACS numbers: 03.70. +k,11.15.-q,12.38.Aw 
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The non-Abelian gauge fields (Yang-Mills fields) [lj], which are introduced by general- 
ization of the local gauge symmetry of the Maxwell electrodynamics from U(l) to SU(n), 
are the basis of the unified theory of electro-weak interactions {U(l) x SU(2)), and strong 
interactions (SU(3)). The quantization of these theories, which is available by now in the lit- 
erature only in the perturbational limit, involves peculiarities and demands an introduction 
of ghost fields to fulfill the unitarity condition 

In the present Letter we will perform the quantization of the Yang-Mills fields by re- 
quiring that the Heisenberg's (quantum mechanical) equations of motion coincide with the 
classical field equations. This approach is closely related to Dirac's method of quantiza- 
tion of a classical system, which consists in appropriate replacement of the Poisson brackets 
by the commutation relations. The Heisenberg's equations of motion and their general- 
ization to the space derivatives (see Eq. (jSJ) for electromagnetic field) can be regarded as 
the differential (local) formulation of the unitarity condition similar to the quantum me- 
chanical systems with a few degrees of freedom. Using the uncertainty relations for electric 
and magnetic field strengthes, which follows from the corresponding commutation relations, 
we estimate the ground state energy of the Yang-Mills fields in the long wave limit as 
&o,k ~ chA^V^ 1 (g 2 /ch) 1 ^ 3 , which depends non analytically on the fine structure constant, 
and on a microscopic cutoff A, which separates the strong and weak coupling limits. In the 
method used in this Letter the nonlinear quantum dynamics of Yang-Mills fields is expressed 
by the simultaneous commutation relations, which allow one to derive using the uncertainty 
relations the expectation values for vacuum fluctuations of field strengthes by fully taking 
into account the non linearities. Our estimate using the uncertainty relations can be consid- 
ered as a simplest version of the Rayleigh-Ritz method of quantum mechanics, and is similar 
from the technical side to the estimate of the ground state energy of interacting oscillators. 
We expect that despite the rather crude derivation the above result is robust and represent 
a lower bound for the vacuum energy of the Yang-Mills fields. 

To elucidate the procedure we will first consider the quantization of the free electromag- 
netic field via Heisenberg's equations of motion 5], and the estimate of the vacuum energy 
using the uncertainty relations. The requirement that the quantum mechanical equations of 
motion for field strengths 
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where 



% = \J rf 3 r(E 2 (r) + B 2 (r)) (2) 

is the Hamilton operator of the free electromagnetic field in Lorentz-Heaviside units, coincide 
with the Maxwell equations 



requires the fulfilment of the following commutation relations for field strengths 

dx' m 



[E l (r',t),E n (r,t)] = 0, 

[B l (r\t),B n (r,t)} = 0. (4) 

The derivatives d/dt in Heisenberg's equations of motion have been replaced by partial 
derivatives d/dt in Maxwell equations. The Heisenberg's equations of motion (CQ) can be 
considered as the differential form of the unitarity condition ensuring the conservation of 
probability during the time evolution of the system. 

The equations of motion with respect to space derivatives are given by 

arm • flnm • 

where 

V = J d 3 rE(v,t) x B(r,t) 

is the momentum of the electromagnetic field. The explicit calculation of the commutators 
in Eqs. ([5]) using the commutation relations (jl]) shows that the quantum mechanical equation 
of motion fl5]) are equivalent to the Maxwell equations 

divE = 0, divB = 0. (6) 

The equations of motion ([5]) can be considered similar to ([1]) as the differential form of the 
unitarity condition ensuring the conservation of the probability by space propagation of the 
electromagnetic field. 

To estimate the vacuum fluctuations of the quantized electromagnetic field using the 
commutation relations (jlj) we follow the consideration for a harmonic oscillator and replace 
p in the energy E = p 2 /2+u 2 q 2 /2 according to p ~ fi/2q, which follows from the commutator 



[q,p] = ih, and results after minimization in the exact expression for the ground state energy 
Eq = hu/2. For the electromagnetic field ones uses the uncertainty relations between the 
strengthes of the electric and magnetic field for a mode propagating in ^-direction in the 



form (AB^) 2 x (Ai±^ fe ) 2 ~ (hk z /2V) 2 . The elimination of the strengthes of the magnetic 



field from the field energy (T5]) and minimization with respect to E^E^ one obtains EkE_k 
= hk/2V, and finally 2 x huk/2 for the vacuum energy, where the factor 2 takes into account 
the number of polarizations. 

The Yang-Mills fields are introduced by requiring the invariance of the Lagrangian with 
respect to the local gauge transformations of the matter field, which in the case of SU(2) 
consists of two components, and transforms according to ip = exp(—ig8 a T a )ip' (in relativistic 
units h = c = 1) with T a = <J a /2, and a a being the Pauli matrices. The local gauge invariance 



demands a coupling to an external Yang-Mi 



the covariant derivative (see for example jf}]) according to 

where the shortcut = W^T a is used. The sum convention over the isotopic index a is 
implied above and in the following. The gauge invariant Lagrangian of the Yang-Mills field 



Is field, which is introduced by replacing <9„ with 



is given by js] 



n }_r^ a CUM 

— ~ ^fiv^a i 



where the components of the tensor of field strengthes are given by 

G% = W - d v W; - ge abc W»Wc. (7) 
The field strengths are expressed through G a in the same way as in electrodynamics 

rpk r~<kO r>k ^klmrylm r~ilm Amk r>k 

Jl/„ = Lr„ , D„ = £ Lr„ , Lx „ = — £ r>„ . 

a a, ' a, ^ a ' a, a. 

The second pair of the classical field equations can be obtained from the Euler-Lagrange 
equations 

dC d dC 



with W£ n = dW^/dx 1 as D v G ilv = 0, and can be written for components £?* as 

d u G^ + ge abc G^W^ = 0. (8) 
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The isotopic indices can be written in equivalent way as subscripts or superscripts. The 
first pair of the field equations (Bianchi identity) follows from the relation between the field 
strengthes and the vector potentials given by Eq. ([70, and reads 

jytQp, + D HQUa + jjvqw = q 

The latter result in the following explicit equations for G^ u 

+ d»G? + d v G°/ + ge ahc G^W a c 
+ge abc GZ°W? + ge ahc G^W v c = 0. (9) 

The energy T 00 and the momentum T ok of the Yang-Mills field are given by js] 

T 00 = \ J d s r (E k E k + B k B k ) , (10) 

l0fc - f d 3 re klm E l a B™ = [ d 3 rGfG ik . (11) 



We now will in quantizing the Yang-Mills fields require that the Heisenberg's equations of 
motion coincide with the classical field equations. This procedure is equivalent to the Dirac's 



method jl| consisting in replacing the Poisson brackets by commutation relations according 
to {A, B} —7- —i/h[A,B]. The field strengths G£ u and vector potentials W£ are operators, 
which generally do not commutate with each others, and are supposed to obey the equations 
(EH). In that case the products G^Wj? should be replaced by (l/2)(G% u Wg + W^G^). 
However, we will not do this explicitly, because our result on the estimate of the energy of 
vacuum fluctuations will not be affected. In order to establish the commutation relations 
for the field strengthes we compare the Heisenberg's equation of motion 

+ yd\>B k (r')[B k (r>),E™(r)} (12) 



with the field equation 



d°G b m0 + d k G b mk - ge bac G a mi W2 = 0. (13) 

The requirement that both equations coincide yields the commutation relations between the 
different isotopic components of the strengths of the electric field 

[E k (r'), E?(r)] = -thge abc 5 km W?5(r' - r), (14) 



and the commutation relation between the strengthes of the electric field and magnetic field 

[B:(r'),E l b (r)] = th5 ab e nlm d' m 5(r> - r) 

- ihge abc e nlm W™5(r' -r). (15) 

The first term on the right-hand side of (IT51) coincides with that for the electromagnetic field 
®- 

To derive the commutation relations for the strengthes of the magnetic field we consider 
the equation of motion 

dB™ 



TtOO nml 

dt ~ 

= \J d 3 r'E k (r>) [E k (r% B?(r) 
+ \J d\'B k b (r>)[B k b (r'),B™(r)] (16) 

and demands that the latter coincides with the field equations (J9]) for a — 0, /J, — i, u — m 

8 

— G a = -a G a - a Gr a - ge abc Lf b w c 

- geabcG^Wl-ge^W™. (17) 



The latter requirement will be fulfilled, if in addition to ( 1151) the commutation relation 
between the strengthes of the magnetic field takes the shape 

[B k (r'), B?(r)] = -ihge abc 6 kn W?6(r' - r). (18) 

The evaluation of the equation of motion 

d m Gf = - % - [T°"\Gt°] (19) 

using the above commutation relation results in 

m G i b ° = -diG[% mi + d m G i b ° - ge abc GfW™ 

- ge abc 5 mt G l a °W< - ge abc GTW° c . (20) 

The part of the latter which is proportional to 5 mi results in the field equation 

diG[ Q + ge^GfW? = 0, (21) 
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whereas the remainder yields the constraint equation 

SabciGf W? - Gl m W° c ) = 0. (22) 
Similarly, the equation of motion 

d m B™ = - % - [T 0m ,B™] , (23) 
decomposes in the field equations 

d m B™ = gs abc B™W™, (24) 
which are associated with terms ~ 8 m% , and the constraint equation 

e ahc {B\W™ - e mn E^W° c ) = 0. (25) 



Therefore, the equations of motion (I19II23P coincide with the corresponding field equations 
only if the constraint equations f)22|25l) are fulfilled. The quantization via equations of 
motion enables one to reduce the complicated quantum dynamics of the Yang-Mills field to 
the commutation relations of the field strengthes (j!4|15|18p at the same time. The constraint 
equations reflect the peculiarities of quantization of Non-Abelian gauge fields. 

In order to check the consistency of the commutation relations fll4lll5|18l) we compare the 
number of commutation relations with the number of vector potentials. The total number of 
commutation relations fll4lll5|ll8l) is equal to 9+4x9+9, while the number of vector potentials 
is 12. The commutation relations between the components of the electric and magnetic field 
and consequently the uncertainties relations of these components are determined according 
to (I14II18P by the zero component of the vector potential W®. Using similar to the case of 
the electromagnetic field |8J] the particular gauge W® = has the consequence that there 
are no uncertainties between different isotopic components of the electric and magnetic 
field strengthes. This choice of gauge decreases the number of constraints. Assuming the 
independence of the expectation values of field strengthes in the ground state on the isotopic 
and the Cartesian indices one decreases further the number of constraints. As a result the 
over-specification of field strengthes will be lifted. 

It is known that the naive canonical quantization of non-Abelian gauge fields is in conflict 



wit 

3, 



i the unitarity condition, which makes necessary an introduction of the ghost fields 
3],{4]. The fixation of gauge to avoid the multiple counting of the field states differing 



only by the gauge occurs by the introduction of Dirac's delta functions under the path 
integral for the vacuum- vacuum transition amplitude, which are the source of the ghost fields. 
The peculiarity of quantization via equations of motion manifests itself in the constraint 
equations (j22|25p . It is beyond the scope of this Letter to discuss the connection between 
both approaches. The constraint equations (I22II25P vanish for expectation values, so that 
they do not affect the treatment of vacuum fluctuations in the following, which is based on 
the uncertainty relations. 

The commutation relation ( !T5|) is the start point to estimate the vacuum fluctuations of 
the energy of the non-Abelian gauge fields, which we will perform here in the long wave limit. 
It is convenient to use for this goal the commutation relations in £;-space. For example, the 
relation between the field strengthes and the vector potentials reads 

G%(k) = ik^w^k) - ik v w;) 

- ge abc J2K( k ') W ^- k - k ')- ( 26 ) 

k' 

We now will derive using the commutation relations the uncertainty relations for field 
strengthes. In considering the limit of small wave vectors we omit the fc-dependent terms 
in commutation relations. As we have mentioned above there are no uncertainties between 
the different isotopic components of field strength, which is the consequence of the choice of 
the gauge W® = 0. The uncertainty relations between the field strengthes of the magnetic 
and electric fields are obtained from the commutation relations (I15D for small k as follows 



Hl{k')El{k) = \e abc \ W 2 c {k + k'), (27) 

and cyclic permutations over the cartesian (upper) indices. The field strength H^(k') in 

/ . \ 1/2 

(I27|) and in the following formulas means H^(k') = lH%(k') 2 ) , and analogous for the 
strengthes of the electric field and vector potentials, where the overbar denotes the quantum 
mechanical average. The elimination of the strengthes of the magnetic field in the energy 
ffTUj) using ( 12"T|) yields after the minimization with respect to the strengthes of the electric 
field 

£o^lhgJ2K( k )> (28) 

k 

where we assumed that the expectation values are independent of the cartesian and isotopic 
indices The factor 9 = 3x3 apply for SU(2). The elimination of the expectation value 
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of the potential in (128]) occurs using fl26|) . and yields 

w -^T^r- (29) 

Note that in the sum over k' in ( 126]) we have approximated the potentials at zero wave 
vectors. Inserting (|29|) into (128]) we arrive finally for the energy £q^ of the fc-mode at 

Further we replace the sum over k by the integral according to Ylk 1 ~ * V I d 3 k/(2ir) 3 = 
VA 3 /6tt 2 , where the microscopic cutoff A in the integral over k has been introduced, and 
obtains in conventional units 

? (3*r 3 ^| (//<*) 1/3 . pi) 

The cutoff A separates the small k from the large ones, where the nonlinearities are expected 
to be less important. The multiplication of ( 13T]) with the factor VA 3 /6n 2 yields the total 
vacuum energy for small k as 

; 2\ 1/3 

47r2 i/3™ y-ft^ 

which is independent of the volume. The numerical prefactor for SU(3) is equal to 3 x 9. 
Note that the energy depends non analytically on the fine structure constant. 

To get a numerical value of the gap energy for SU(2) we use a w = g 2 /ch ~ 0.004, 
set A equal to the inverse of the interaction radius of weak interactions A = m^c 2 /he with 
mw — 80 GeV/ c 2 being the mass of the H^-boson, and obtain £su(2) — 7.2 GeV. The estimate 
for SU(3) with a s ~ 1 and setting 1/A equal to the radius of the proton ~ 0.84 x 10~ 15 m 
yields £su(3) — 0.95 GeV, which is of order of magnitude of the nucleon mass. 

To conclude, we have derived using the quantization of Yang-Mills fields via the quantum 
mechanical equations of motion the non linear commutation relations for the field strengthes 
of non-Abelian gauge fields at the same time. The equations of motion with respect to the 
space coordinates ( fl9lT23l coincide with the corresponding classical field equations only by 
imposing the constraints ( )22l25l) . From the commutation relations between the magnetic 
and electric fields we derive the uncertainty relations for expectation values of vacuum 
fluctuations of electric and magnetic field strengthes, and use the latter to estimate the 



£o-JwW^) . (82! 



energy of vacuum fluctuations of the Non-Abelian gauge fields in the long wave limit. The 
vacuum energy does not vanish for k — > i.e. it possesses a gap, which is expected to 
represent a lower bound of the vacuum energy of the Yang-Mills fields. 



[1] C. N. Yang and R. Mills, Phys. Rev. 96, 191 (1954). 
[2] R.P. Feynman, Acta Physica Polonica, 24, 697 (1963). 
[3] L.D. Faddeev and V.N. Popov, Phys. Lett. B 25, 29 (1967). 
[4] B. De Witt, Phys. Rev. 162, 1195, 1239 (1967). 

[5] P. A. M. Dirac, Directions in Physics, John Wiley & Sons, New York, 1978. 

[6] C. Ytzykson and J.B. Zuber, Quantum Field Theory, Vol. 2, MacGraw Hill, New York, 1980. 

[7] P. A. M. Dirac, The Principles of Quantum Mechanics, Clarendon Press, Oxford, 1958. 

[8] V.B. Berestetskii, E.M. Lifshitz, L.P. Pitaevskii, Quantum Electrodynamics. Vol. 4 (2nd ed.), 

Butterworth-Heinemann, 1982. 
[9] Note that in performing the minimization we neglect the relation between the field strengthes 

and the vector potential. 



10 



